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Abstract 

In this paper we discuss limits on magnetic fields that could have 
been present at nucleosynthesis time. We considered several effects 
that could be relevant modifing light elements relic abundances. They 
include: changes in reaction rates, mass shifts due to strong and elec- 
tromagnetic interactions, variation of the expansion rate of the Uni- 
verse due to both the magnetic field energy density and the increasing 
of the electrons density in overcritical magnetic fields. We find that 
the latter is the main effect. It was not taken into account in previous 
calculations. The allowed field intensity at the end of nucleosynthesis 
(T = 1 x 10 9 °K) is B < 3 x 10 10 Gauss. 
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1 Introduction 

Amongst the many uncertainties in the Early Universe environment, the pos- 
sibility that large, constant magnetic fields, existed over macroscopic scales, 
is a fascinating possibility. Since the Early Universe is believed to be a perfect 
conductor, magnetic lines get thinned out by the ratio 

Bx = (BW 

The presence of primordial fields over galactic scales, when extrapolated back 
can give, under different assumptions, very large fields indeed. 
A typical present day galactic field B = 10~ 6 G, can grow, scaled as dictated 
by Eq.(P alone to be as large as 10 14 Gauss. These extrapolations are very 
doubtful since dynamo effects may have enlarged significantly present fields. 
Nevertheless, knowing the allowed fields at a given epoch and limiting its 
value at another one can give important dynamical restrictions. 
Recently Vachaspati has shown that large magnetic fields, of the order 
of B = (mw) 2 = 10 24 Gauss, should be generated at the electroweak phase 
transition time. His argument seems quite general. It only depends on the 
finiteness of the horizon and embedding electromagnetism in a larger field. 
The size of a patch at electroweak time, is of order \ = N{rnw)~ l - Evolving 
the size of the patch is model dependent. More precisely, the coefficient N in 
[|I|] is a function of the scale. If the scale is determined by today's size then 
N = 10 13 . This argument predicts a field of about a Gauss at nucleosynthesis 
time. 

The size of the patch at the nucleosynthesis time might allow for larger fields, 
though we have no real reliable model for the field evolution. 
In this paper we address ourselves to the allowed magnetic fields at nucle- 
osynthesis time, without discussing their origin. We give a detailed analysis 
of the influence of the fields on the main quantities that can act to modify the 
relic elements abundance ratios: reaction rates, masses of the participants, 
electron energy densities and magnetic field energy density. We then obtain, 
given the present errors in these relative abundances, the upper limits the 
fields can take in regions large compared to the reaction scale but possibly 
much smaller than the horizon at that time. In the first section we discuss 
briefly the impact of the magnetic field on these elements of the calculation. 
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In the next we describe the standard nucleosynthesis calculations in this light. 
In the final section we discuss the constraints and their origin and compare 
with existing calculations. 



2 Weak reaction rates in the presence of mag- 
netic fields 



The main weak processes which act to determine the n/p ratio during the 
the primordial nucleosynthesis are 
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The rate for two body scattering reactions in a medium may be written in 
the form 
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where Pi is the four momentum, Ei is the energy and f\ is the number density 
of each particle species. All processes in Eq. (a,b and c) have the same 
amplitude 
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where a = gA/gv — — —1-262. Without any external magnetic field the total 
rate of the processes that convert neutrons to protons is 
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where - = G ^"^3 ^ and q, e and (fi e are respectively the neutron-proton 
mass difference, the electron energy and the electron chemical potential all 
expressed, in units of m e . We assume the neutrinos chemical potential to be 
vanishing. 

The total rate for the p — > n processes can be obtained changing the sign 
of q in Eq.(||). 

An external magnetic field take us to modify Eq. (|j) due to the following 
effects. 

a) The dispersion relation of charged particles propagating trough a mag- 
netic field is modified with respect to the free-field case. In fact, their 4- 
momentum is in this case p = p(B = 0) + qA , where q is the charge of the par- 
ticle and the vector potential A(r) is related to the field by A(r) = |r x B. 
Assuming B along the z axis, the expressions for energies of electrons, pro- 
tons and neutrons are respectively 
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In the above, n denotes the Landau level, s = ±1 indicates whether the spin 
is along or opposed to the field direction, and g p = 5.58 and g n = —3.82 are 
the Lande g-factors. The QED correction to the electron energy, k, has been 
first computed by Schwinger P|. For magnetic fields larger than ~ 10 13 G 
this correction is 
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For smaller field intensity k has negligible effects on our calculations and 
we disregarded it. The effects of the field on the QCD ground state have 
been parametrized via a field dependent nucleon mass [0] , as we are going to 
discuss below. 

Neither the neutron or the neutrino have quantized levels, though the 
neutron has an electromagnetic interaction energy. The neutrino is totally 
inert vis a vis electromagnetism. 
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b) The number of available states for a particle obeying Eq. (4) becomes, 
for every value of n and s, is|| 

This changes the phase space of the processes we are interested to consider. 

c) Since the occupation number and the energy of states with opposite 
spin projections is not the same in a magnetic field, the spin sum of the 
square amplitude needs to be weighted by the appropriate spin dependent 
Fermi distributions. 

Nucleosynthesis take place in a range of temperatures 0.1 < T < 10 MeV, 
hence nucleons are nonrelativistic. Therefore nucleon distribution functions 
are given by 

f N (s = ±1) = (1 + e^)- 1 (10) 

where 

Since during the nucleosynthesis rriN ^> T, and momenta are also small 
compared to nucleon mass, can be safely approximated by 1/2. This is 
not the case for electrons. In this case we have 

fe(s) = (l + e E ^ T y 1 (12) 

where the relativistic expression for the electron energy Eq.([|) is used. 

As a consequence, the integral for the leptonic momentum space in neu- 
tron (3 decay is modified to 

_L £ r % ['■■<" % eB _ m _ fe) (13) 

27T n=0 J-oo 1E V J—p e ,z(n) 2E e 

where N c is largest integer n such that p e , z {n) 2 = Q 2 — m 2 e — 2neB is positive 
and Q 2 = m 2 n — m 2 

d) The nucleon masses are affected by very strong magnetic fields. The 
change in effective phase space is 

A = 0.12fi N B-M n + M p + f(B) . (14) 
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The function f(B) gives the rate of mass change due to colour forces being 
affected by the field. For nucleons || the main change is the chiral condensate 
growth, which because of the different quark content of protons and neutrons 
makes the proton mass grow faster 0. Though the sign is certain, vacuum 
pairs of zero helicity get more bound in the presence of a B field, the size 
of the effect is model dependent. We have calculated its influence using the 
weakest and strongest reasonable field dependence and we find the effect 
always small for fields below 10 18 Gauss. 

Having established that hadronic mass changes will not affect nucleosyn- 
thesis we drop these effects from the equations altogether. 

Taking into account the remaining effects we computed the total rate for 
the weak processes converting neutrons to protons in an external magnetic 
field. The result is 
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where 7 = B/B c and B c 
critical magnetic field. 



m 2 e /e = 4.4 x 10 13 G is usually defined to be the 



Equation fll5"|) coincides with the result of Matese and O'Connell |5| and 
Cheng et al. || in the limit in which the QED correction k goes to zero. 
Although the quantitative effects of this term on the nucleosynthesis pre- 
dictions are subdominant, we stress that disregarding it when the field is 
overcritical it leads to an unstable electron ground state, thus to unphysical 
results. 

The main effect of the magnetic field is due to the modification of the 
electron phase space. Eq.(|i~5|) is correct in the weak field limit, when 7 1 
and T n ^ p (B) reduces to Eq.(f|) in the B = limit. 

In Fig.(l) we present the dependence of T n ^ p as a function of temperature 
for some values of 7. 

For large values of 7 and fixed temperature, the total rate grows like 
7. Increasing the temperature the relevant contribution to the integrals in 
Eq.(|T5|) comes from the high energy part of the electron spectrum. Since 
the limit e — > 00 is equivalent to the limit 7 — >• in Eq.(O) this explains 
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why the ratio T n ^ p (B)/T n ^ p (0) goes to one when T 3> m e . Although 
the global rate of the inverse process T p ^ n also increases with B, it remains 
suppressed by a factor exp(—Q(B)/T) with respect to T n ^ p . Thus the effect 
of a strong magnetic field would be to reduce the final number of neutrons in 
the Universe, i.e. the relic 4 He abundance, if only the correction to the weak 
rates is taken into account. 



3 The effects of B on the expansion rate 

Owing to exponential dependence of the (n/p) equilibrium ratio on the tem- 
perature, the relic relative abundances of light elements depends crucially 
on the freeze-out T F temperature of the weak processes that keep protons 
and neutrons in chemical equilibrium |7|. This temperature is essentially 
determined by the condition 

IW7» = H(T F ) (16) 

where H is the expansion rate of the Universe. 

It is evident that besides the rate of the weak processes we need to pay 
attention to the effects of the magnetic field on H. If no cosmological constant 
is present, the expansion rate is determined by the Einstein equation 

H\T) = ^p(T) . (17) 

where p(T) is the total energy density of the Universe. In the case that no 
magnetic field is present p(T) is given by the sum of the energy density of 
all the particle species in thermal equilibrium with the primordial plasma 

p{T) = p 7 (T) + Pe (T) + p u {T) + p b (T) (18) 

where the subscripts 7, e, u, b stand, respectively, for photons, electrons, 
the three species of neutrinos and baryons, including their respective anti- 
particles. In our case, since the magnetic field has energy density Pb{T) = 
£>(T) 2 /87r this term also needs to be added to Eq. flT8|) . Since we have mag- 
netic flux conservation in the plasma 

B oc R- 2 oc T 2 
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the energy density of the magnetic field has the same temperature dependence 
as the energy density of the radiation. 

This new contribution to p(T) will dominate over the other terms in 
Eq.flD if 

B(T = 10 llo K) >10 16 G. (19) 

We assumed the pressure associated to the random magnetic field to be 
zero in average. Although a novanishing mean pressure is also possible for 
random magnetic fields M, our final conclusions are not affected also taking 
this pressure into account. 

The presence of a nonvanishing pb is not the only effect that modifies 
the expansion rate of the Universe. The energy density of charged particles 
in the primordial plasma is also affected. In the previous section we have 
showed how the electron dispersion relation and the electron phase-space are 
modified by the magnetic field. Using Eqs.(|) and (|9]) we get the electron 
energy density in function of 7 

eB 00 r 
Pe(T) = fj £ / dp z E e (s)f e = 

Z7F n=0,s J 

— >(2 - S n0 ) / dee met . 

2tt 2 ^0 V v /i+2(n+i) 7 +« _ K )2 _ 1 _ 2 (n + 1)7 1 + e~ +** 

(20) 

The reader can see from Fig. (2) that the effect of an overcritical magnetic 
given by Eq.(J19|) is to increase the electron energy density roughly linearly 
with the field intensity. The same of Eq . (|20|) is valid for positrons once 
the sign of the chemical potential <p e has been changed. In analogy with 
what happen for the reaction rates, this effect becomes less important at 
high temperatures if B is left fixed. However it is a very relevant effect when 
T < 1 MeV. Although the field intensity, hence the correction to p e , decreases 

like T 2 , we are going to show that this is the main effect on the primordial 
nucleosynthesis predictions. 

4 Results and conclusions 

In the previous sections we have shown that the existence of large magnetic 
fields during the primordial nucleosynthesis affects the final light elements 
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Table I. Predictions of light element abundances at the end of the primordial 
nucleosynthesis are given for several values of the magnetic field intensity 
given at the temperature T = 10 11 °K. 
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relative abundances via two main effects: 

a) the increasing of the weak reaction rates and 

b) the increasing of the expansion rate of the Universe. 

These are competing effects. In fact, whereas the former tend to reduce the 
(n/p) freeze-out temperature, hence the final abundance of 4 He, the latter 
acts in the opposite direction. 

We modified the standard nucleosynthesis code |J to take into account 
all the relevant effects, as well as other effects that eventually we neglected 
as irrelevant. 

Since our aim is to get an upper limit to the magnetic field intensity, we 
adjusted the value of the baryon photon ratio rj in order to get the minimal 
4 He relic abundance prediction compatible with observations [10] in the free- 
field case. In Tab. I we present our predictions for some light elements relic 
abundances. As it is evident at a glance, limits on magnetic fields are totally 
controlled by the 4 He abundance. 

Other elements reach forbidden values only at very high fields, in that 
case the effect of mass changes due to colour forces will also be important. 

The increase of the 4 He relic abundance with the field intensity reveals 
that the effect of B on the expansion rate is the most relevant. Regarding this 
point we agree with the qualitative conclusion of Matese and O'Connell 



8 



and disagree with the opposite conclusion of Cheng et al. [II]]. Mainly, we do 
not understand how they can reconcile the claim that the dominant effects of 
the magnetic field are those arising from modification of the reaction rates, 
with the growing of the relic 4 He that they get increasing B. 

Furthermore, in both refs. |J and the effect of the magnetic field on 



the electron and positron energy density was not considered. Leaving only 
this effect on in our code, we checked that this is indeed the most relevant. 

We showed that this is indeed the main effect as long as the field intensity 
at the beginning of the nucleosynthesis is smaller than 10 16 G. 

Since the observational upper bound for the 4 He relic abundance is Y p < 
0.245 we conclude that the average intensity of a random magnetic field at 
the temperature of T = 1 x 10 11 °K (beginning of nucleosynthesis) must 
be less than 3 x 10 14 G or, equivalently, B{T = 10 9o K) < 3 x 10 10 G 
(end of nucleosynthesis). Vachaspati |TJ predicts a magnetic field strength of 
~ 10 11 G, on the smallest coherence region of the field (N = 1), at the end 
of nucleosynthesis. This extreme assumption (N = 1) is ruled out by our 
limits. 

Assuming the field continue to rescale according to Eq.(l) (perhaps not a 
reasonable assumption), our results imply that the intergalactic field is less 
than ~ 3 x 1CT 7 G at present. 

The other light elements relic abundances are less affected than He by 
the magnetic field and we do not use them to get constraints. However, it 
is interesting to observe the behaviour of Deuterium and 3 He abundances 
versus the initial magnetic field. Although they increase with the field at 
the beginning of nucleosynthesis their relic abundances follow the opposite 
behaviour. This can be understood since the rates of the processes converting 
Deuterium and 3 He to 4 He are proportional to the initial abundances Yd or 
F3 He 0- The greater are the rates smaller are the freeze-out temperatures 
for these light elements. Thus we expect smaller D and 3 He relic abundances 
even if the relic 4 He increases. 
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gure Caption 

1. The neutron-depletion rate T n ^ p , normalized to the free- field rate, is 
plotted as a function of the temperature for some constant values of 7. 

2. Electron energy densities are displyed as a functions of 7. The lower 
curve corresponds to the temperature value T = 0.5 MeV; the middle 
one to T = 1 MeV and the upper one to T = 5 MeV. 

3. Electron and magnetic field energy densities temperature dependence 
are compared for two different values of the initial magnetic field. 
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